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iD ' Abstract 



We suggest a rigorous generalization of the pseudoscalar mesons mixing description in 
SU (3) basis. It is shown that the appearance of extra massive state nicknamed in the paper 
as glueball is unavoidable in any scheme with more than one angle. In this framework we 
develop the dispersive approach to Abelian axial anomaly of isoscalar non-singlet current. 
Combining it with the analysis of experimental data of charmonium radiative decays ratio we 
get the number of quite strict constraints for mixing parameters. Our analysis favors the equal 
values of axial currents coupling constants which may be considered as a manifestation of 
SU (3) symmetry and possible violation of chiral symmetry based predictions. 
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1 Introduction 

The problem of pseudoscalar meson states mixing has been under intense theoretical and experi- 
mental study for many years. This topic attracts a lot of interest due to its close relation to such 
fundamental phenomena as quantum anomaly, chiral symmetry breaking and study of exotic states 
like glueball, which can elucidate the essential features of QCD. 

During the last decade, a large amount of new experimental data on mesons has been collected. 
New data from the forthcoming experiments at COMPASS, BES-III, GlueX (JLAB) and at the 
upgraded facilities FAIR (GSI) and MAMI might allow a complete quantitative verification of the 
different mixing schemes and approaches. 

A number of analyzes of mixing in i] — rf system based on different processes have been 
performed in the last decades and the mixing angles in the range —10° -. — 24° were obtained (see, 
for example (TJ - IH ). The analysis of the axial anomaly generated decays 77(7/) — > 27 was also 
performed in [|9l (in the framework of ChPT) and ifTOl . and the estimation 9 = —20° -r 25° was 
obtained. 



*e-mail: klopot@theor.jinr.ru 
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■'■e-mail: teryaev@theor.jinr.ru 
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The approach with one mixing angle dominated for decades in the studies of the 77-77' system. 
However, in the recent years there was a rise of interest in the mixing schemes with two and three 
angles. The theoretical ground of this was based on the observation that taking into account the 
chiral anomaly through perturbative expansion in ChPT can lead to the introduction of two mixing 
angles in the description of the 77-77' system [fTTL IT2T | . Some analyzes of various decay processes 
were also performed in this scheme lTT3l IT4l IT5l IT6l ITTTl . 

The theoretical analysis of the mixing in the pseudoscalar sector base either on SU(3) or quark 
basis. The last was introduced by T. Feldmann and P. Kroll [13] and was widely used in the last 
years. 

For us it happens to be more convenient to construct and use the rigorous generalization of 
SU(3) basis similar to the mixing of massive neutrinos. This is because we use the dispersive 
approach to axial anomaly [fT8l (see also [19] for the review) to find some model-independent and 
precise restriction on the mixing parameters (angles and decay constants) in the different schemes 
(with one and more angles). The combination of our approach with certain processes leads to quite 
strict predictions. We consider this paper as a first step and application of our approach to the 
whole set of processes is to be done. 

This paper is organized in the following way. In Section 2 we introduce the general mixing 
scheme in SU (3) basis. More specifically, we derive the expression for (non-square) matrix of 
coupling constants generalizing the approach offered by B. Ioffe and M. Shifman [|20l |2T |. In 
Section 3 we consider the dispersive approach to axial anomaly for the isoscalar non-singlet ax- 
ial current J^ 5 taking into account the possible contributions of higher mass states. In this way 
we derive our main equation while in Section 4 we supplement it by the analysis of the ratio of 
radiative decays of J/ty. These two constraints happen to be sufficient to provide strong bounds 
for mixing parameters. In Section 5 we consider the reduction of the general scheme with three 
angles to some currently popular particular (two-angle) schemes Il22l - [|28l . Finally, in Section 6 
we summarize the results of numerical analysis and discuss the possible implications for SU (3) 
and chiral symmetries. 



2 Mixing 



We start with a vector of physical pseudoscalar fields consisting of the fields of lightest pseu- 
doscalar mesons and other fields: 



V 

rf 
G 

w 



(1) 



We do not need to specify the physical nature of other components with higher masses, the lower 
of which G is probably a glueball. Let us also introduce, following [|20ll2TI . a set of SU(3) fields 
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<fk{k = 3, 8, 0) and other (singlet) fields gf. 
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(2) 



and corresponding states \Pk) and \Pk>)- The three upper fields ip 3 , (p 8 , y? diagonalize the matrix 
elements of axial currents J^ 5 = q^^Ti^^q : 



(0\J l jPk) = iSikfkQn , I, fe = 3,8,0. 



All other states are orthogonal to these currents: 



<o|jyjv> = o. 



(3) 



(4) 



At the same time all the corresponding fields enter the mass term in the effective Lagrangian with 
a generally speaking non-diagonal mass matrix M : 



AC = -<P T M<P 

2 

This formula immediately implies the generalized PCAC relation: 



(5) 



d^ut = FM<£, (6) 

where 
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(7) 

\0 f ... 0/ 

F is a matrix of decay constants^ defined in ©. 

In order to proceed from initial SU (3) fields <P to physical mass fields <P the unitary matrix U 
is introduced 



<P = \J<P (8) 

which diagonalizes the mass matrix 

UMU T = M = diag(mlo,ml, m^,m%, . . .), (9) 

where m^, m,,, m^/ and are the masses of the it, r\, rj mesons and glueall state G respectively. 

'Note, that matrix of decay constants F is non-square expressing the fact that generally the number of SU(3) 
currents is less then the number of all possible states involved in mixing. The similar situation takes place (see e.g. 
ll29ll ) in one of the extensions of the Standard Model — neutrino mixing scenario involving sterile neutrinos. 
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Simple transformations of Eq.© reads: 





cV m5 = FU T M<£ (10) 

This formula is close to those obtained in EOl l2"T]| (in the limit of small mixing). When the decay 
constants are equal this is reduced to formula (3.40) in 11301 . 

Taking into account the well-known smallness of n° mixing with 77, 77' sector EOl |2T1 [36 1 and 
neglecting all higher contributions we restrict our consideration to three physical states 77, 77', G and 
two currents J^ 5 , J° 5 . Then the divergencies of the axial currents (recall, that G is a first mass state 
heavier than 77'): 

£8 = (o fa oKff ^ ° ' m 

Exploring the mentioned similarity of meson and lepton mixing, we use a well-known general 
Euler parametrization for the mixing matri)@ U: 

/ C 8 C 3 — C0S3S8 — C 3 S 8 — CgCQSs S3S0 
1 I »:>, <■-.', I I r,, -.s (l I I S 8 C 8 I = S 3 C 8 + C 3 C S 8 -s 3 s 8 + C 3 C 8 C -c 3 s 

\ s 8 s c 8 s c 

(12) 

As soon as in the chiral limit J^ 5 should be conserved, Eq.dTTb obviously implies, that the 
coefficients in front of the terms m^, rriL, m 2 G should decrease at least as (m, / m v i : g) 2 ■ 

The generic matrix (fl2|) can be reduced to different particular cases. The pure 77 — 77' mixing 
with no other admixtures corresponds to the case #0 = in Eq. (fT2~)) . This is the so-called one- 
angle mixing scheme with the mixing angle 9 = 9 8 + 6> 3 . As it was mentioned above, that this 
scheme is not sufficient for description of the full set of experiments. 

One can easily see from (fTTI) . (fT2l) that the schemes with more than one mixing parameter 
unavoidably require introduction of new states for the mixing matrix to be unitary. Indeed, from 
these equations one can see that in general case we have 3 different angles and in some particular 
cases one can reduce the number of angles by one and to get the schemes with two different angles 
(put 9 3 = or 9 8 = 0). Subsequently the schemes with 3 and 2 angles will be considered in details 
in sections 4 and 5. 

In this work we sequentially use SU(3) basis. For our purposes this basis is more preferable 
since in the next Section we will consider non-singlet isoscalar axial current J^ 5 which is free from 
non-Abelian anomaly. The transition to quark basis <3> q = ((«« + dd)/y/2, ss,g) which is also 
widely used in literature can be performed by multiplying by additional rotation matrix V: 

-v/273 0\ 

$ = V$ q , V = VT73 . (13) 

V 1/ 



2 We use notation c, = cosOi , s, = sinOi 
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3 Dispersive approach to axial anomaly 



In our paper the dispersive form of the anomaly sum rule will be extensively used, so we remind 
briefly the main points of this approach (see e.g. review lfT9ll for details). 

Consider a matrix element of a transition of the axial current to two photons with momenta p 
and p' 

T^(p,p') = ( P ,p'\J^\0) . (14) 

The general form of T^ a g for a case p 2 = p' 2 can be represented in terms of structure functions 
(form factors): 



T^ a /3(p,p') = F 1 (q 2 )q fM e a(3p(T ppp' CT + -F2(q 2 )[-^e^ p(T p p p' a -— e mpa p p p' a - e^^p-p 1 )^ , (15) 

L p p 

where q = p + p'. The functions Fi(q 2 ), F 2 (q 2 ) can be described by dispersion relations with no 
subtractions and anomaly condition in QCD results in a sum rule: 

oo 

J ImF 1 (q 2 )dq 2 = 2aN c Y,e 2 q , (16) 
o 

where e q are quark electric charges and iV c is the number of colors. This sum rule was proved 
I13T1 for the general case p 2 ^ p' 2 and earlier in Il32l and 11331 for the cases p 2 < 0, m = and 
p 2 = p' 2 respectively. Notice that in QCD this equation does not have any perturbative corrections 
11341 , and it is expected that it does not have any non-perturbative corrections also due to 't Hooft's 
consistency principle 11351 . It will be important for us that at q 2 — > oo the function ImFi(q 2 ) 
decreases as 1/g 4 . Note also that the relation (fT6l) contains only mass -independent terms, which is 
especially important for the 8th component of the axial current J^ 5 containing strange quarks: 

jjS = ^(«7m75W + hulsd - 257^755) (17) 

The general sum rule (fl~6l) takes the form: 

00 

J Im F l {q 2 )dq 2 = -j=a(e 2 u + e 2 d - 2e 2 )N c = ^-a , (18) 


where e u = 2/3, = e s = —1/3, N c = 3. 

Consider now a particular case of pure r] — r]' mixing, where no other mixing states are taken 
into account. Recall that this case has been studied for a long time, various approaches were 
considered and the numbers for the mixing angles in the range — (10 4- 24) were obtained. The 
approach based on the dispersive representation of axial anomaly was introduced for 77 — 7r° mixing 
in 11361 and used in [|371 where r] — r]' mixing was considered in assumption of small mixing angle). 

In order to separate the form factor F 1 (q 2 ), multiply T Ma( g(p, p') by q p /q 2 . Then taking the 
imaginary part of Fi(q 2 ), using the expression for d p J^ 5 from Eq. (fTTI) and saturating the matrix 
element with the r], 77' states we get: 
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Figure 1 : Mixing angle Q x as a function of the decay constant / 8 in the one-angle mixing scheme. 
Dashed curves correspond to the uncertainties of the experimental data input. Horizontal dot- 
dashed line indicates the / 8 = 1.28/^ level. 



ImFi(q 2 

h 



Im q^\ 



(27 



t(8) 



0) 



(2 7 | [m^( Cl )+m^'(si)] I 0) 



wf 8 [A v 5(q 2 - m?)(ci) + A v ,5(q 2 - ml)( 8l )] 



(19) 



Note, that if we have included higher resonances to this equation, they are expected to be sup- 
pressed as \jm\ es by virtue of the mentioned above asymptotical behavior of F x {q 2 ) oc l/q\ This 
(approximate) independence on higher resonances together with (exact) quark mass independence 
of the anomaly relation (fT6l) may be an indication of some connection between these two effects. 

If we employ the sum rule (fTHT) , we obtain a simple equation: 



ci + /3si = f, 



where 




a 2 m^ 



967r 3 r r 



f2 ' 
-»2 7 J 8 



H\a 2 

64vr r ' 



(20) 



(21) 



(22) 



For numerical evaluation of the 77 — 7/ mixing angle in the latter case, put / 8 = 1.28/^ , 
/ w = 130.4 MeV, m v = 547.85 MeV, m' r] = 957.78 MeV, T^ 27 = 0.51 keV, V n ,^ 2l = 4.30 keV. 
The value for / 8 is taken from the chiral perturbation theory calculations, other numbers from PDG 
Review 2008 [f38l . The mixing angle, responsible for the mixing of the 77 — 7/ system appears to be 
0i = —22. 1° ±1.5°. The dependence of the mixing angle on / 8 is shown on the FigJT] Remarkably, 
the anomaly sum rule fixes the mixing angle (provided we know / 8 ) in this case. 

From our point of view one of the advantages of our approach is a high accuracy due to a high 
accuracy of the anomalous sum rule (TTSl for octet axial current. Let us stress that at this stage 
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we avoid consideration of the anomaly relation for the singlet axial current which contains the 
contributions of gluons, direct instantons and topological effects (see e.g. 11391 , [[40 J ) . As a result 
it appears to be unnecessary to include f to our analysis. 

The current result moderately agrees with our earlier analysis in the small mixing angle approx- 
imation [37]. It is not completely trivial that our result is also in a good agreement with previous 
analysis done in pioneering papers fl2j[5]|. This is because Eq.(T20l) happens to follow also from 
the (non-dispersive) anomaly equations used in [5j . The similar result of dispersive and non- 
dispersive (local) approaches is in some sense natural taking into account that we omit the higher 
contributions with controlled accuracy 0(l/m 4 ). At the same time the mentioned approaches 
actually use the anomalous divergency for singlet current which we do not need. 

Now let us consider a generic case involving glueball admixtures. Performing the same opera- 
tions for the r] — rf — G system we get: 

ImF 1 (q 2 ) = Imq fl ±{2 1 \J$\0) = 

- — (27 | [m 2 r](c 8 c 3 - c s 3 s 8 ) + m 2 ,r]'(s 3 c 8 + c 3 c s 8 ) + m|Gs 8 s Q ] I °) = 

tt/s [A v 5(q 2 - m 2 ) (c 8 c 3 - c s 3 s 8 ) + A v ,5(q 2 - m 2 , ) (s 3 c 8 + c 3 c s 8 ) + A G 5 (q 2 - m 2 G )(s 8 s )} (23) 
The final equation following from the anomaly sum rule (TT8l) for the 77 — rf — G system is: 

(c 8 c 3 - c s 3 s 8 ) + p(s 3 c 8 + c 3 c s 8 ) + 7(s 8 s ) = f , (24) 

where 



i = 

Let us summarize the situation with a theoretical accuracy of the Eq. (1241) . As we have pointed 
out, the anomaly sum rule (fT8b has no a s -corrections^. The possible contributions from higher 
states may come as the additional terms in the l.h.s. of Eq. (|24l) . As it was discussed before, the 
asymptotic behavior Fi(q 2 ) is proportional to 1/g 4 at large q 2 and a sort of quark-hadron duality 
implies that higher resonances should be suppressed as (m^//m res ) 4 . 
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a 2 m^ 1 
96tt 3 IV 27 if 



G^2j m r i . 



r 



3 ' 



nr 



—La 2 
64tt v 



(25) 



(26) 



4 J/^ radiative decay ratio 

Eq. (l20l) provide an exact constraint but contains too many free parameters. As an additional 
experimental constraint we, following 0, J21 use the data of decay ratio Rj/qr = (T(J/^f) — > 

rfi)/(r{j/*)^m). 

3 However, this is not the case for the singlet axial current J° 5 . 
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As it was pointed out in [HTTl the radiative decays J/ty — > 77(^)7 ^ dominated by non- 
perturbative gluonic matrix elements, and the ratio of the decay rates Rj/^ = (T(J/^) — > 
r} lr ))/{Y{J/^!) — > 777) can be expressed as follows: 



(0\GG\ i) 



(0 I GG I 77) 



(27) 



where = M J /y(l—rn% l , T)l JMj,y)/2. The advantage of such a choice is an expected smallness 
of perturbative and non-perturbative corrections. 

The divergencies of singlet and octet components of axial current in terms of quark fields can 
be written as: 



9^*5 = -^=(m u M7 M 7 5 u + rnddj^d - 2m s s7 / /y 5 s), (28) 



1 1 3ct 

d^J% = -j={rn u u^f^u + m d d-f^ 5 d + m s s^f^ 5 s) + -^-j=-^-GG (29) 

Following [0, neglect the contribution of u- and d- quark masses, then the matrix elements of 
the anomaly term between the vacuum an 77, vj states are : 



^(0 1 gg 1 V ) = (o 1 1 V ) + ±(0 1 a„jg> 1 V ), oo) 



^(0 I GG I V ') = (0 I dpJ® I rf) + -L(0 I d^jjS I t/> (31) 
Using Eq. (Qj), (ED), ©, (ED we deduce: 

/o(-s 3 s 8 + c 3 c 8 c ) + ^/ 8 (s 3 c 8 + c 3 c s 8 ) 



R 



™* /o(-c 3 s 8 - c 8 c s 3 ) + ^/ 8 (c 8 c 3 - c s 3 s 8 ) 



(32) 



Let us note that in obtaining this equation only operator relations for anomalies (1281 [291 were 
used and one did not need to express the 77(77') mesons fields in the terms of divergencies of singlet 
(and octet) axial currents. 

If we use Eq.(|32l) for the case of one-angle mixing scheme (6 = 0), for usual choice / 8 = 
1.28/t,., /o = 1.1 /•„■ and corresponding angle value obtained from anomalous dispersive relation 
(|20l) [Ox = 9 3 + 9 S = —22.1°) we find Rj/q, = 2.2 which is in serious discrepancy with the 
experimental value Rj/q, = 4.8 ± 0.6 [|38l . Substituting the experimental value of R into (|32|) and 
using anomaly equation (1201) one can get the dependencies / 8 (#i), /o(#i) (Fig- HI). 

As soon as we accept, that fs>fo> fir (f° r different kind of justification see e.g. [fT4l . IfTTTl ). 
it follows from Fig. [2l that the only possibility is f 8 ~ f ~ 1 (#1 ~ —18°), which is quite far 
from the chiral perturbation theory expectations 0. Taking into account all experimental errors (the 
dominant contribution being provided by that of Rj/q>) one get Fig. [3] where the effects of these 
errors are indicated by shaded areas. From this figure it is clear that the maximal allowed value 



4 These values appear to be preferable also in other mixing schemes as we will see later. 
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e u [degrees] #i [degrees] 



Figure 2: Relationship between decay constants 
fs, fo and mixing angle 9\ in the one-angle mix- 
ing scheme 



Figure 3: Relationship between decay constants 
fs, fo and mixing angle 9\ in the one-angle mix- 
ing scheme, rescaled. Dashed curves denote ex- 
perimental uncertainties, meshed area indicates 
the range of parameters satisfying the condition 

/8 > fo > /tt 



of fs is fs = 1.2f n (= fo). Let us note that these values correspond to minimal allowed value of 
Rj/q, = 4.2. Here the importance of more accurate experimental value of Rj/q, is already clear. 

Let us pass to the more elaborated schemes with more than one mixing angle which were of- 
fered in [TT1IT31 . In the most general scheme with 3 angles introduced in Section 2 it is convenient 
to rewrite d32]) in terms of 9 1 = 9$ + 9 3 , 9 2 = 9 8 - 9 3 : 



m\, fo(c 1 -c 2 + Co(c 1 + c 2 )) + ^f 8 ( y s 1 -s 2 + c ( y s 1 + s 2 )) 
J * ™% /o(-«i - s 2 - c (si - s 2 )) + Tj/stci + c 2 - c (c 2 - d)) 

The angle 9± has an explicit physical sense. From the definition of the mixing matrix U (fT2l) 
it is obvious that the angle 9% describes the overlap in the r] — v( system with accuracy ~ 9^/2 
and coincides with their mixing angle at 9 — > 0. It is reasonable to suppose that the glueball 
contribution to the anomaly sum rule for non-singlet current (|24l) is rather small. This doesn't 
necessarily mean the extreme smallness of rj — G mixing angle itself but rather the cumulative 
effect of smallness of mixing angles and large mass m G . That's why we neglect the last term in 
the l.h.s. of (|24l) and rewrite it in terms of 9±, 9 2 : 

(ci + c 2 - c (c 2 - ci)) + /3(si -s 2 + c (si + s 2 )) = 2f . (34) 

The solutions of (|33l (upper curves, red online), and (|34l ) (lower curves, blue online) are shown 
on Fig. |4]for customary choice of decay constants fs = 1.28f w , fo = 1.1/*- and three different 
angles 20°, 32°, 40°. We see from these figures that at 9 = 20° there is no intersection. The first 
common solution appears at larger angle 9 = 32° when the curves touch each other. As 9 grows, 
at 9 > 32° two different solutions appear. While for these two solutions 9 2 are significantly 
different, the solutions for 9\ are limited to relatively narrow region 9\ = — 13° -. — 18°. This 
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^^^^/s/ fir 
fo/ fn 


1.28 


1.1 


1.0 


1.28 


24^ 






1.1 


qo+a.t) 
OZ -4.5 


17.5^ 2 




1.0 


35^ 


23.5_ 6 5 


10+ 9 

iu -10 



Table 1: Tree-angle mixing scheme. Minimal possible mixing angle #0 (in degrees) for different 
values of decay constants f and / 8 

is not surprising because of the physical sense of 9i mentioned above 0. Taking into account the 
experimental uncertainties we see that the minimal allowed value of 9 slightly decreases to 27°. 

All this numerical analysis showed that for / 8 = 1.28, f = 1.1 we need a substantial glueball 
admixture 9 > 27°. The minimal value of 9 decreases only if f 8 decreases. The results of analysis 
for different f 8 and f are presented in Table [H from which one can make an important conclusion 
that relatively small glueball admixture even within experimental uncertainties is possible only for 
f 8 = fo (and most probably f 8 = f = f w .) 




-150-100 -50 50 100 150 -150-100 -50 50 100 150 -150-100 -50 50 100 150 

e 2 [degrees] e 2 [degrees] 6 2 [degrees] 



(a) O = 20° (b) O = 32° (c) O = 40° 

Figure 4: General (tree-angle) mixing scheme. Solutions of Rj/y equation (|33l) (upper. red online 
curves) and anomaly condition (|34l) (lower, blue online curves) for / 8 = 1.28/^, fo = 1.1 f w and 
different choice of 9 . The dashed curves indicate experimental uncertainties. 



5 Two-angle mixing schemes 

Up to a moment our analysis was quite general. It is instructive to consider some particular cases 
which are currently discussed in literature. Starting from our general scheme, one can easily 
perform a reduction to partial mixing schemes with 2 angles. Here we have 2 distinct cases when 
ind 

i) # 3 = or 

ii) £ 8 = 0. 

5 Other solutions not shown on the figures appear only at large O > 70° which clearly has no physical sense. 
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1. The case 6 3 = 0. 

This case was introduced in 112211 . and was widely used by KLOE collaboration in a set of recent 
papers Il23ll - ll25l . This choice clearly means that there is no mixing of rj with additional scalar 
state noted in © as g, i.e. r\ is a mixture of 4> 8 and </>o only. In this case our equations (|24)) . (|32l) 
are simplified as follows: 

c 8 + f3c s 8 + <ys s 8 = £, (35) 



m\, fo(c 8 c ) + -^fsicoss) 



m\ / (-s 8 ) + ^/s(c 8 ) 



(36) 



Note that 8 in this scheme defines r\ — rj mixing, as one can easily see from Eq. (12). In 
further analysis we suppose that 7 cannot exceed 2 for any reasonable values of Y G ^ 2 ^ (i-e. 
T G ^2-y/ m G IS 4Tri-^2-f/ m ^)- This restriction corresponds to the assumption that 2-photon de- 
cay widthes of pseudoscalar mesons grow like the third power of their masses, or in other words 
the glueball coupling to quarks to be of the same order as for the meson octet states 

The results of numerical analysis are shown on Fig. [5l These figures show the dependence of 
the glueball contribution 7 on the angle 9 8 for different values of f 8 and f . 

On Fig. |5Ta)fl the dependence y(6 8 ) is shown for / 8 = 1.28/^ and f = (1.1, 1.28)/^. The 



dotted lines corresponds to experimental uncertainties, the uncertainty for R being dominant. From 
this figure one can see that for any fo<f 8 — 1-28/^ we get 7 > 2. One can achieve 7 ~ 2 only 
for f ~ 1.28 ~ / 8 (at the lower value of Rj/^ = 4.2 ). The corresponding values of mixing 
angles are 9 8 = -(14 4 17)°, 9 = (12 4 30)° . 

On the Fig. |5(b)| the case / 8 = 1.1 for two choices of fo = (1.1, 1.0) is shown. Again, the 
reasonable values 7 < 2 is achievable only for f ~ 1.1 ~ / 8 . The corresponding values of 
mixing parameters are 8 8 = —(12 -i- 18)°, 9 = (5 4- 35)° 



And, at last, consider the choice fo = f 8 = 1-0/V shown on Fig. |5(c)| The corresponding 
values of mixing parameters are 9 8 = — (10 -j- 18)°, 6 = (5 4- 37)° Note that in this case the region 
of relatively small 7 < 1 and 9 ~ 5 are achievable. 

Let us notice that for / 8 > 1.28/^ the minimal value of 7 is growing, e.g. for large N c value 
/ 8 = 1.34/,,- (and any f ), 7 > 2 within experimental errors of Rj/y. 

We can conclude, that in this scheme the glueball admixture is bounded from below to minimal 
value 6> > 5° for any / 8 > f > f w . The relatively small glueball admixture 9 < 10° is possible 
only for / 8 ~ / - U- 

Let us stress, that here (like in general 3-angle mixing scheme) the values / 8 ~ f is much 
more preferable. All these results directly follow from the specific choice of mixing scheme, axial 
anomaly condition and Rj/y. 

It is instructive to compare our results with those obtained in 11251 , where the analysis of 
BR{<p 7/7) / BR{4> — * 777) together with several vector meson radiative decays to pseudoscalars 



6 Note that in these papers authors use quark basis, see (TOT l, The connection between notations of mixing angles in 
ll25ll and our notations: 8q — —4>g, #8 = 4>p — a, where a = arctan a/2/3 ~ 54.7°. 

7 If it grows with mass according to the naive dimensional arguments, the glueball width will be even smaller. We 
will discuss the relation between coupling constants and decay widthes later. 

8 We limit ourselves to negative 9$ since this mixing angle in this scheme clearly have a sense of r\ — rf mixing 
angle. 
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was performed in the same mixing scheme. 

If we substitute their values of 9 8 = (— 14.5 ± 1), 9 = 7.5° to our dispersive relation and take 
into account possible sources of inaccuracy we can see that we directly get the dependence 7(/ 8 ) 
shown on Fig. [6] The large slope is due to small mixing angles 9 and 9 8 . 




ft B [degrees] a s [degrees] e [degreesl 



(a) / 8 = 1.28/, (b) / 8 = 1.1/, (c) fa = 1-0/. 



Figure 5: Glueball contribution parameter 7 as a function of mixing angle 9 8 for different decay 
constants. Dashed curves indicate errors arising from experimental data. 

One can see that the reasonable value 7 < 2 corresponds to f$ ~ 1((1.00 — 1.05)/^). Thus 
anomaly constraint contradicts the preferred value [llj f s = 1.34/^ (/ 8 = 1.28/^) which is 
used in [25 J. However, the constant f s enters only the expression for radiative decay width 
7/ -> 2 7 (Eq.1.4 in [|25l). The change of f s from f s = 1.34/ 7r (/ 8 = 1.28/,) to the value 
fs = fn(fs = fir) in this equation may be compensated by the change of within the claimed 
accuracy. Therefore at the current level of accuracy the analysis of KLOE is compatible with 
ours. The combination of the KLOE results with both our constrains (i.e. anomaly and Rj/y, (|35l 
l36l) leads to conclusion that fs = fo = fn which is in agreement with the analysis presented above. 

2. The case 9 8 = 0. 

One can consider another particular case of general mixing supposing 9 8 = in (fT2l) . This 
kind of a two-angle mixing scheme was used recently in 11281 . This scheme implies no glueball 
coupling to J^ 5 , which means that the glueball contribution to anomaly relation is exactly zero. 
Therefore, the anomaly relation is exactly the same as for the one-angle mixing scheme so the 
dependence f 8 (9 3 ) is the same and shown on Fig. \T\ Recall, that this relation is precise due to 
reasons, mentioned in Section 3. 

The relation for J/\& radiative decay ratio for this case can be directly obtained from (1321 with 
#8 = 0. This relation together with axial anomaly relation allow us to get the dependence of the 
ratio of decay constants fo/fs as a function of 9 (see Fig. [7]). Stripes bounded by dashed 
lines show the effects of experimental uncertainties. From Fig. |7]it is clear that the solution at 
/s = 1.28/, contradicts to the condition /o//s < 1 discussed above for any #0. The only solution 
compatible with /o//s < 1 is the one at fs = (1 — 1.05)/, (the left stripe), the corresponding 
6*o = 0° 4- 14°, while the angle #3 (having the sense of the 77 — 7/ overlap) can vary in a very narrow 
region 9 3 = -17° 4 -18°. 

So one can conclude that this scheme even within current experimental uncertainties demands 
/s = /o = fir an d small glueball admixture 9 < 14°. 

9 As it was noted in Section 4 the ratio / / fs should not exceed 1 . 
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Figure 6: Two-angle mixing scheme, case 6 3 = 
0. Glueball contribution parameter 7 as a func- 
tion of decay constant / 8 (in units of f n ) for 
6$ = (-14.5 ± 1)°, O = 7.5° 
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Figure 7: Two-angle mixing scheme, case #3 = 
0. The dependence of the mixing angle Oq on 
the ratio of decay constants fo / /g together with 
experimental uncertainties. 



Let us now examine the numerical results of [28J. Implementation of our procedure with their 
angle <p = 42.4° (6 3 = -12.3) leads to / 8 = 0.8 and contradicts to f 8 /f n > 1.0 being the robust 
expectation from the ChPT. 

As soon as the accuracy of the extracted [|28l parameters is not available at the moment it is 
difficult to conclude whether this disagreement is statistically significant. 



6 Summary 

In this paper we presented the detailed analysis of the role of Abelian axial anomaly in the mix- 
ing of both light and heavy pseudoscalar states. We found that the anomaly imposes the severe 
constraints for the meson couplings and mixing angles. There is also a delicate interplay between 
pseudoscalar state, which we call glueball without specifying its nature, and light mesons. 

We offered a (new, to our best knowledge) rigorous approach in SU(3) basis to the consider- 
ation of mixing of pseudoscalar mesons. Our approach is quite similar to that for the mixing of 
massive neutrino, where the number of states with a definite mass may exceed the number of flavor 
states. In this sense the appearance of extra singlet mesons is similar to the role of sterile neutrinos. 

We use the dispersive representation of axial anomaly for J^ 5 . The advantage of such rep- 
resentation is the high and controlled accuracy due to suppression of higher resonances and its 
independence on the quark masses. The use of SU(3) basis allows us to limit ourself to more 
theoretically clear case of non- singlet current receiving contribution only from Abelian anomaly. 
Let us note, that our equation (|20l) may be formally derived in the other pioneering approaches 
using standard local anomaly relations (see e.g. B21 in the one-angle mixing scheme. 

As a supplementary input we use experimental data for the ratio of radiative decays Rj/y = 
(r(J/^) — > 7/7) / (r( J/ty) — > 777) which is theoretically safe and provides additional restriction. 

We found that for the one-angle mixing scheme the only reasonable solution is / 8 — /o = 
f-K, @i = —18°. 
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We proved that any scheme with more than one angle unavoidably demands additional singlet 
admixture and considered the general scheme with 3 angles with one lowest additional singlet 
state, which we denoted as a glueball without any specification of its nature. 

It was found, that the combination of our two inputs (anomaly relation and Rj^ data) leads to 
a rather strict constraints. The values / 8 = 1.28f n inevitably leads to a large glueball admixture 
9 ~ 30°. At the same time, the reasonably small glueball admixture 9 < 10° is possible only for 
/o ~ fs = l-O/n-, which is far from expectations, based on chiral perturbation theory. 

We checked this general observation by considering some particular cases, when the general 
mixing matrix (with 3 angles) is reduced to the i] — rj' — G mixing scheme with two angles. 

The first one is the mixing scheme with 6*3 = which was investigated in a set of recent 
papers by KLOE collaboration. We concluded that application of our approach to this scheme 
unavoidably leads to fs/ fo ~ 1 for any reasonable values of glueball two-photon decay width. 
Moreover, combining our constraints with the angles 9 8 = 14.5°, 9 = 7.5° obtained in KLOE 
analysis based on the additional set of decays, we again found that f 8 rs f = f n . These values are 
still consistent with the results of their fits within their accuracy. 

We also consider another two-angle mixing scheme where 9 8 = 0. In this case the only possible 
solution compatible with constraint /§ > / is fs = fo = fn with high accuracy. 

As a result, in all the considered schemes the relation / 8 ~ / ( ~ most likely) holds. This 
marks a sort of new manifestation of SU(3) symmetry and at the same time the possible violation 
of chiral perturbation theory expectations. The possible origin of such a symmetry pattern may be 
the smallness of the strange quark mass (squared) with respect to nucleon one while it is still much 
larger than (genuine) higher twist parameters and may be treated sometimes as a heavy one [42J. 

The significant progress may be achieved by the more accurate determination of Rj/q,, in par- 
ticular, at BES-III accelerator, complementing its vast program (see [|43il, Section 17.2.2). While 
our conclusion / 8 ~ f is more robust and is valid for all values of Rj/y within current experi- 
mental limits, the stronger result / 8 = /„. may be questioned by the more accurate data. Therefore 
these measurements will provide a new test of chiral perturbation theory predictions. 

Although our approach already provided the important new constraint for the analysis in pseu- 
doscalar channel, the global fit exploiting the dispersive representation of axial anomaly remains 
to be done. 
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